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ABSTRACT 

The brightest galaxy in a dark matter halo is expected to reside at rest at the center of 
the halo. In this paper we test this 'Central Galaxy Paradigm' using group catalogues 
extracted from the Two-Degree Field Galaxy Redshift Survey (2dFGRS) and the 
Sloan Digital Sky Survey (SDSS). For each group we compute a parameter TZ, which 
is defined as the difference between the velocity of the brightest group galaxy and 
the average velocity of the other group members (hereafter satellites) , normalized by 
the unbiased estimator of the velocity dispersion of the satellite galaxies. Since the 
redshift surveys suffer from incompleteness effects, and the group selection criterion 
unavoidably selects interlopers, a proper comparison between data and model needs 
to take this into account. To this extent we use detailed mock galaxy redshift surveys, 
which are analyzed in exactly the same way as the data, thus allowing for a fair 
comparison. We show that the central galaxy paradigm is inconsistent with the data 
at high confidence, and that instead the brightest halo galaxies have a specific kinetic 
energy that is about 25 percent of that of the satellites. This indicates that either 
central galaxies reside at the minimum of the dark matter potential, but that the halo 
itself is not yet fully relaxed, or, that the halo is relaxed, but that the central galaxy 
oscillates in its potential well. The former is consistent with the fact that we find a 
weak hint that the velocity bias of brightest halo galaxies is larger in more massive 
haloes, while the latter may be indicative of cored, rather than cusped, dark matter 
haloes. We discuss several implications of these findings, including mass estimates 
based on satellite kinematics, strong gravitational lensing, halo occupation models, 
and the frequency and longevity of lopsidedness in disk galaxies. 

Key words: galaxies: halos — galaxies: kinematics and dynamics — dark matter — 
methods: statistical 



1 INTRODUCTION 

In the standard picture of galaxy formation, hot gas in viri- 
alized dark matter haloes cools and accumulates at the cen- 
ter of the potential well, where it forms a galaxy (White & 
Rees 1978). During the hierarchical build up of larger and 
larger structures, haloes with their 'central' galaxies are ac- 
cumulated by even larger haloes. At that point the halo be- 
comes a subhalo, and the central galaxy becomes a satellite 
galaxy. In the standard picture, it is envisioned that a satel- 
lite galaxy no longer accretes hot gas, which instead is only 
accreted by the galaxy in the center of the potential well 
(e.g., Kauffmann, White & Guiderdoni 1993; Somerville & 
Primack 1999; Cole et al. 2000). Since this central galaxy 
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therefore continues to grow, it is expected to be the bright- 
est, most massive galaxy in a halo. This is further assured 
by the fact that any other massive galaxy would quickly sink 
to the center of the potential well by dynamical friction to 
merge with the central galaxy, thus producing an even more 
massive central galaxy. Therefore, according to the standard 
paradigm, the brightest galaxy in a halo will reside at rest 
at the center of the potential well. Note that this is clearly a 
statistical statement, as it does not necessarily hold for each 
individual system (e.g., non-virialized, strongly interacting 
systems). Hereafter, we will refer to this paradigm as the 
'Central Galaxy Paradigm' (CGP), and use the terms 'cen- 
tral galaxy' and 'brightest halo galaxy' without distinction. 

The CGP plays an important role in various areas of as- 
trophysics. For example, attempts to measure halo masses 
from the kinematics of satellite galaxies, are always based 
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on the general assumption that the 'host' galaxy is located 
at rest at the center of a relaxed halo (e.g., Zaritsky et 
al. 1993, McKay et al. 2002; Brainerd & Specian 2003; Prada 
et al. 2003; van den Bosch et al. 2004). This assumption 
is also used in virtually all mass models of strong gravita- 
tional lenses. On the other hand, the observed frequency 
and longevity of lopsidedness in disk galaxies (e.g., Richter 
& Sancisi 1994; Zaritsky & Rix 1997) is often interpreted as 
evidence for an actual offset between galaxy and halo (e.g., 
Levine & Sparke 1998). The central galaxy paradigm also 
plays a role in halo occupation modeling, where assump- 
tions have to be made regarding the spatial distribution of 
galaxies in haloes in order to compute the galaxy-galaxy 
correlation function on small scales (e.g., Scoccimarro et 
al. 2001; Berlind & Weinberg 2002; Yang, Mo & van den 
Bosch 2003; van den Bosch, Yang & Mo 2003; Maglioc- 
chetti & Porciani 2003; Tinker et al. 2004; Zehavi et al. 2004; 
Zheng et al. 2004). A statistic that is particularly sensitive 
to whether the brightest halo galaxies reside at the center 
or not is the cross correlation between dark matter haloes 
and galaxies (see Yang et al. 2005, in preparation). 

This special dynamical status of the brightest galaxy in 
a halo has been tested for the special class of cD galaxies. 
Jones et al. (1979) have shown that cDs are located at the 
peak of the cluster X-ray emission, while Quintana & Lawrie 
(1982) used the kinematics of cluster galaxies to argue that 
cDs are at rest with respect to the cluster. Although this 
is in agreement with the CGP, more recent studies have re- 
vealed various cases in which the cD galaxy has a significant 
peculiar velocity with respect to the mean velocity of the 
other cluster members (e.g., Sharpies, Ellis & Gray 1988; 
Hill et al. 1988; Zabludoff, Huchra & Geller 1990; Oegerle 
& Hill 1994, 2001). Applying a similar study to a dozen 
poor groups, Zabludoff & Mulchaey (1998) and Muchaey 
& Zabludoff (1998) found that the position of the bright- 
est galaxy in each group is indistinguishable from that of 
the group center or from the center of the X-ray emission. 
To our knowledge, however, the central galaxy paradigm 
has never been tested for a statistically significant sample 
of dark matter haloes that span a wide range in masses. In 
this paper we use data from the Two-Degree Field Galaxy 
Redshift Survey (2dFGRS, Colless et al. 2001) and the Sloan 
Digital Sky Survey (SDSS; York et al. 2000) to directly test 
whether the brightest galaxies in dark matter haloes are lo- 
cated at rest at the center of their potential well. We show 
that, although the brightest halo galaxies are clearly segre- 
gated with respect to the other galaxies in the same halo, 
they have a typical specific kinetic energy that is about ~ 25 
percent of that of the satellite galaxies, and that the CGP 
is ruled out at a high level of confidence. 

This paper is organized as follows. In Section 2 we 
present a statistic that can be used to test the CGP, which 
we apply to the 2dFGRS and SDSS in Section 3. In Section 4 
we describe a simple model for the velocity and spatial bias 
of the brightest halo galaxies, which we use in Section 5 to 
construct detailed mock galaxy redshift surveys of the 2dF- 
GRS. In Section 6 we compare these mocks with the data in 
order to constrain the phase-space parameters of brightest 
halo galaxies. Section 7 discusses various implications of our 
results, and we summarize our conclusions in Section 8. 



2 DYNAMICAL SIGNATURE OF CENTRAL 
GALAXIES 

Observationally, the only kinematic information that is 
available to test the central galaxy paradigm are the line- 
of-sight velocities obtained from redshifts. In what follows 
we use v c to refer to the line-of-sight velocity of the bright- 
est halo galaxy, and Vi is the line-of-sight velocity of the 
i th satellite galaxy. In addition we define the difference 
AV — Vs — v c between the mean velocity of the satellite 
galaxies (v s — jj- Vi ) an d that of the central galaxy. 

If the CGP is correct and Vi follows a Gaussian distribution 
with velocity dispersion a s , the probability that a halo with 
N s satellite galaxies has a value of AV is given by 



P{AV)dAV 



2%o 



exp 



(AVl 
2a 2 



dAV, 



(1) 



with a — a s / ' \/ N s 
the parameter 



R 



Therefore, in principle, one could define 



(2) 



and test the CGP by checking whether R follows a normal 
distribution with zero mean and unit variance. However, the 
velocity dispersion a s is generally unknown, and we have to 
use its unbiased estimator 



(3) 



instead. This allows us to define the modified parameter 



(4) 



If the null-hypothesis of the CGP is correct, 1Z should follow 
a Student t-distribution with v = N s — 1 degrees of freedom. 
Note that P V (TZ) approaches a normal distribution with zero 
mean and unit variance in the limit N s — > oo. 

The applicability of this 'TJ-test' is strongly related to 
the ability to find those galaxies that belong to the same 
dark matter halo. To this extent we use the halo-based 
galaxy group finder developed by Yang et al. (2005a), which 
has been optimized for this task. Although this group finder 
is well tested and calibrated, it is not perfect. In particular, 
because of redshift errors and redshift space distortions, it 
is unavoidable that one selects interlopers (galaxies that are 
not associated with the same halo). The expectation value 
of \v s — v c \ will be larger for an interloper than for a true 
satellite. As long as the interloper is fainter than the bright- 
est galaxy in the group (halo) to which it is assigned, its 
impact on 1Z may be small, as it affects both the numerator 
and the denominator. However, if the interloper is brighter 
than all true group members, \Ti\ will typically be severely 
overestimated. Another problem is related to the fact that 
the 2dFGRS and SDSS suffer from various incompleteness 
effects. If the actual brightest halo galaxy is missed (i.e., is 
not present in the survey) , 1Z will be measured with respect 
to a satellite galaxy, which again will bias \1Z\ high. The 
presence of interlopers and incompleteness effects, therefore, 
tend to create excessive wings in the 1Z distribution. A com- 
parison with the Student t-distribution might then give the 
wrong impression that the null-hypothesis is rejected. Since 
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Figure 1. A comparison of the cumulative distributions of \TZ\ of 2dFGRS (red, solid lines) and SDSS (blue, dashed lines) groups. 
Results are shown for four intervals in a s , indicated in square brackets in each panel. The KS probability, Pks, that both distributions 
are drawn from the same distribution is also indicated. Note that the P(\1Z\) from 2dFGRS and SDSS are in excellent agreement with 
each other. 



the typical occupation numbers of haloes are small, this ef- 
fect can be very strong, as we demonstrate in Section 6. To 
circumvent these problems, we compare the ^-distributions 
obtained from groups in the 2dFGRS and SDSS against 
those obtained from groups extracted from detailed mock 
galaxy redshift surveys, which suffer from interlopers and 
incompleteness effects to the same extent as the real data. 

3 APPLICATION TO THE 2DFGRS AND SDSS 
3.1 Group Selection 

The 7?.-test described above requires a selection of galax- 
ies that belong to the same dark matter halo. In Yang 
et al. (2005a, hereafter YMBJ) we developed a halo-based 



galaxy group finder, that is optimized for this task. Here we 
give a brief description of this group finder, and refer the 
interested reader to YMBJ for details. 

The basic idea behind our group finder is similar to that 
of the matched filter algorithm developed by Postman et 
al. (1996), although it also makes use of the galaxy kinemat- 
ics. The group finder starts with an assumed mass-to-light 
ratio to assign a tentative mass to each potential group, iden- 
tified using the friends-of-friends (FOF) method. This mass 
is used to estimate the size and velocity dispersion of the un- 
derlying halo that hosts the group, which in turn is used to 
determine group membership (in redshift space). This pro- 
cedure is iterated until no further changes occur in group 
memberships. Using detailed mock galaxy redshift surveys, 
the performance of our group finder has been tested in terms 
of completeness of true members and contamination by in- 
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Figure 2. The velocity bias (left-hand panel) and spatial bias (right-hand panel) of central galaxies as function of the parameter / CC n, 
which expresses the characteristic scale of the radial distribution of central galaxies in terms of the characteristic scale of the NFW 
density distribution (see Section 4). Results are shown for three values of the halo concentration parameter c, as indicated. 



terlopers. The average completeness of individual groups is 
~ 90 percent and with only ~ 20 percent interlopers. Fur- 
thermore, the resulting group catalogue is insensitive to the 
initial assumption regarding the mass-to-light ratios, and is 
more successful than the conventional FOF method in as- 
sociating galaxies according to their common dark matter 
haloes. 

3.2 The 2dFGRS 

We use the final, public data release from the 2dFGRS, re- 
stricting ourselves only to galaxies with redshifts 0.01 < z < 
0.20 in the North Galactic Pole and South Galactic Pole 
subsamples with a redshift quality parameter q > 3 and a 
redshift completeness c > 0.8. This leaves a grand total of 
151,280 galaxies with a sky coverage of ~ 1125 deg 2 . The 
typical rms redshift and magnitude errors are 85kms _1 and 
0.15 mag, respectively (Colless et al. 2001). Absolute mag- 
nitudes for galaxies in the 2dFGRS are computed using the 
K-corrections of Madgwick et al. (2002). 

Application of the halo-based group finder to this 
galaxy sample, yields a group catalogue consisting of 77, 708 
systems. Detailed information regarding the clustering prop- 
erties and galaxy occupation statistics of these groups can 
be found in Yang et al. (2005a,b,c). In what follows we re- 
strict our analyzes to the 2502 groups in this catalogue with 
four members or more. 

3.3 The SDSS 

In addition to the 2dFGRS, we also use data from the SDSS. 
In particular, we use the New York University Value-Added 
Galaxy Catalogue (NYU-VAGC) t, described in detail in 

t http://wassup.physics.nyu.edU/vagc/#download 



Blanton et al. (2004). The NYU-VAGC is based on the SDSS 
Data Release 2 (Abazajian et al. 2004), but with an inde- 
pendent set of significantly improved reductions. From this 
catalogue we select all galaxies in the Main Galaxy Sample, 
which has an extinction corrected Petrosian magnitude limit 
of r = 18. We prune this sample to those galaxies in the red- 
shift range 0.01 < z < 0.20 and with a redshift completeness 
c > 0.7. This leaves a grand total of 184,425 galaxies with 
a sky coverage of ~ 1950 deg 2 . From this SDSS sample, we 
construct a group catalogue that contains 102, 935 systems. 
A more detailed description of this catalogue will be pre- 
sented in Weinmann et al. (2005, in preparation) . As for the 
2dFGRS, we restrict our analysis to the groups with four 
members or more, of which there are 3473 in our catalogue. 

3.4 Comparison of 2dFGRS with SDSS 

For each group in both the 2dFGRS and SDSS catalogues 
described above we compute 1Z. Fig. 1 plots the cumula- 
tive distributions of \1Z\ for both surveys. In the upper-left 
panel we plot the distributions using all groups in the range 
50 kms -1 < a s < 1000 kms -1 , with a s the unbiased es- 
timator of g s (equation [3]). In the other three panels we 
plot P(< \1Z\) for three sub-samples (the values in square 
brackets indicate the range in a s used, in kms -1 ). Over- 
all the agreement between SDSS and 2dFGRS is extremely 
good. To make the comparison more quantitative, we use 
the Kolmogorov-Smirnov (hereafter KS) test to compute 
the probability Pks that both P(|72.|) are drawn from the 
same distribution. The resulting probabilities are indicated 
in each panel. These confirm what can already be inferred 
by eye, namely that both 7£-distributions are consistent with 
each other. Given this good agreement between both data 
sets, we only concentrate on the 2dFGRS in what follows. 
The main reason for choosing this survey over the SDSS is 
that we have accurate mocks for the 2dFGRS that have been 
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Figure 3. The cumulative distributions of \1Z\ obtained from MGRSs Afo.o (left-hand panel) and Mi.o (right-hand panel), in which the 
brightest halo galaxies have a velocity bias of 6 vc i = and b ve i = 1, respectively. Solid, dashed and dotted curves correspond to group 
samples with 50kms — 1 < a s < 200 kms -1 , 200 kms -1 < a s < 350 kms -1 , and 350 kms -1 < a a < 1000 kms -1 , respectively. The 
gray area indicates the area bounded by Student distributions with 3 and 9 degrees of freedom. In the ideal case without interlopers, 
the P(< |7£|) of Mo.o should fall in this range. The fact that they don't illustrates the impact of interlopers and completeness effects, 
and emphasizes the importance of using MGRSs for a fair comparison with the data. Finally, the fact that the P(< \TZ\) of Afo.o and 
Mi.o are significantly different illustrates that the 7^-test does have the ability to constrain the phase-space parameters of brightest halo 
galaxies. 



well tested. Given the good agreement between 2dFGRS and 
SDSS, we argue that any result based on the former will also 
hold for the latter. 



o3»(r) = VX-£s ( - ) ( 1+ - ) Z{r/r„) 



7(c) Vr 

with f(x) = ln(l + x) — x/(l + x) and 



(6) 



4 MODELING VELOCITY BIAS OF CENTRAL 
GALAXIES 

The main goal of this paper is to use the 1Z distributions 
presented above in order to constrain the phase space pa- 
rameters of brightest halo galaxies. We will express these in 
terms of their spatial and velocity bias with respect to the 
satellites. If the null- hypothesis of the CGP is correct, both 
the spatial and the velocity bias should equal zero. In order 
to model these biases, and to incorporate them in the mock 
redshift surveys that we will use for comparison with the 
data, we proceed as follows. 

We assume that each dark matter halo has an NFW 
(Navarro, Frenk & White 1997) density distribution, Pdm(r), 
with virial radius r v i r , characteristic scale radius r s , and con- 
centration parameter c = r v i v /r s . Assuming haloes to be 
spherical and isotropic, the local, one-dimensional velocity 
dispersion follows from solving the Jeans equation 



2 I 
°"dm( 



Pdm(r) 



Pdx 



(r )-Q^i r ) dr 



(•») 



with ^(r) the gravitational potential (Binney & Tremaine 
1987). Using that d^/dr = GM(r)/r 2 and defining the 
virial velocity Kdr = \J GM/r v i T we obtain 



Ay) 



t'HI + tY 



(7) 



The halo-averaged velocity dispersion is given by 
4tv 



Pdm(r) a d m(r) r dr 



M 



/ 3 (c) J (1 + y) 



(8) 



(cf. van den Bosch et al. 2004). 

Throughout this paper, we assume that the iV sa t satel- 
lite galaxies in a halo of mass M follow a number den- 
sity distribution n sa t(r) = (iV sa t / 'M) /5 dm (r) , i.e., there is 
no spatial bias between satellite galaxies and dark matter 
particles. As shown in van den Bosch et al. (2005), this is 
consistent with the observed radial distribution of satellite 
galaxies in the 2dFGRS. If we further assume that the satel- 
lites are in isotropic equilibrium, it also follows that there 
is no velocity bias between the satellites and the dark mat- 
ter, neither globally [i.e. (<7 S at)M = ((Tdm)jif] nor locally [i.e. 
ffsat(r) = cr dm (r)]. 

When stacking all haloes of a given mass, we assume 



© 2000 RAS, MNRAS 000, 1-13 



6 F.C. van den Bosch et al. 




2 4 6 0246 

|R| [Rl 




0246 0246 



|R| |R| 

Figure 4. The cumulative distributions of \1Z\ obtained from the groups in the 2dFGRS (gray dots), compared with those obtained from 
three of our MGRSs, which only differ in their value of b ve i, as indicated in the upper left panel. Results are shown for four intervals in 
o"s , indicated in square brackets in each panel. The number of 2dF groups in each of the three subsamples is indicated. 



that their brightest halo galaxies follow a number density 
distribution given by a Hernquist profile^: 



la 1 
p ce n{r) oc — -- 



2n r (r + a) 3 



(9) 



This implies a probability distribution for r of the central 
galaxies of 



P cen (r)dr = 2 



(r + a) 



rdr 



(10) 



In order to parameterize the characteristic radius a in terms 
of that of the dark matter halo, we define the parameter 
/con = a/r s . A brightest halo galaxy at a halo-centric radius 



* The choice for this particular distribution is not motivated by 
any physical considerations, other than the fact that it is well 
behaved, both at r = and at r — > oo 



r has an isotropic velocity dispersion 



2 / \ I 



„ -oSo' 



V 2 
v,r /( 



with 

J(y) 



/(r)dr 



T 3 (/cen+T)3 



This implies a halo-averaged velocity dispersion of 
J rvlr Pcon(r) er ccn (r) r 2 dr 



(0"cen)j' 



J rv " Pcen(r) r 2 dr 



(12) 



" d f(c) fcen J (Un + y)^ dy ' (13) 

which allows us to define the velocity bias of brightest halo 
galaxies as 6 ve i = {a ccn } / {a d m} = {a ce n} / (a Ba , t ) ■ In addition 
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to the velocity bias, we define the spatial bias as fe rad = 
(r C on}/(rdm) = (r- cen )/(r S at), where the expectation value for 
the radius follows from 



Jo™ P(r)r 3 dr 
Jo V ' r P( r ) r2 dr 



(14) 



For an NFW density distribution with concentration c, this 
reduces to 



(r-dm) = 



(2 + c)/(l + c)-(2/c)ln(l + c) 
ln(l + c) - c/(l + c) 



(15) 



or (rdm) = 0.41r v i r for c = 10. 

Fig. 2 plots 6 ve i (left-hand panel) and & rat j (right-hand 
panel) as function of / ccn for three values of the halo concen- 
tration parameter c. In the limit / con — > 0, the probability 
distribution P CO n(7") becomes a Dirac delta function. This 
implies that the central galaxy is sitting still at the center of 
the dark matter halo (i.e., the null- hypothesis of the CGP), 
so that 6 ve i = fcrad = 0. Increasing / ccn increases the proba- 
bility to find the brightest halo galaxy at larger halo-centric 
radii, which corresponds to a larger velocity bias. Note, how- 
ever, that b vc i never approaches unity, which is due to the 
fact that pcen(r) can not be made to match pdm(r) for any 
value of /ccn- Typically b vo i 3> fe ra d, which is a reflection of 
the 'depth' of the NFW potential. For example, for a ve- 
locity bias of & V ei = 0.5 (i.e., corresponding to a specific 
kinetic energy that is one quarter of that of the satellites) 
the radial bias is 6 rad ~ 0.07 (assuming c = 10). Combining 
this with (15) implies an expectation value for the offset 
of the central galaxy from the dark matter distribution of 
(r ccn ) ~ 0.03r v i r . For a Milky- Way sized system this corre- 
sponds to ~ 5 kpc, comparable to the characteristic radius 
(scale length) of the galaxy itself. 

In what follows, we construct a set of mock galaxy red- 
shift surveys (hereafter MGRSs) for different values of b ve \, 
and compare the ^-distributions of their groups against 
those of the 2dFGRS, which are statistically identical to 
those of the SDSS, in an attempt to constrain 6 ve i. 



5 MOCK GALAXY REDSHIFT SURVEYS 

We construct MGRSs by populating dark matter haloes with 
galaxies of different luminosities. The distribution of dark 
matter haloes is obtained from a set of large A-body simu- 
lations (dark matter only) for a ACDM 'concordance' cos- 
mology with Q, m = 0.3, = 0.7, h = 0.7 and a$ = 0.9. 
In this paper we use two simulations with N = 512 3 parti- 
cles each, which are described in more detail in Jing & Suto 
(2002). The simulations have periodic boundary conditions 
and box sizes of Lbox = 100/i -1 Mpc (hereafter Lioo) and 
ibox = 300/i _1 Mpc (hereafter Z/30o). We follow Yang et 
al. (2004) and replicate the L 30 o box on a 4 x 4 x 4 grid. The 
central 2x2x2 boxes, are replaced by a stack of 6 x 6 x 6 
Lioo boxes, and the virtual observer is placed at the center 
(see Fig. 11 in Yang et al. 2004). This stacking geometry cir- 
cumvents incompleteness problems in the mock survey due 
to insufficient mass resolution of the L300 simulations, and 
allows us to reach the desired depth of z max = 0.20 in all 
directions. 

Dark matter haloes are identified using the standard 
FOF algorithm with a linking length of 0.2 times the mean 



inter-particle separation. Unbound haloes and haloes with 
less than 10 particles are removed from the sample. In Yang 
et al. (2004) we have shown that the resulting halo mass 
functions are in excellent agreement with the analytical halo 
mass function of Sheth, Mo & Tormen (2001). 



5.1 Populating Haloes with Galaxies 

In order to populate the dark matter haloes with galaxies 
of different luminosities, we use the conditional luminosity 
function (hereafter CLF), $(L|M), which gives the average 
number of galaxies of luminosity L that resides in a halo of 
mass M. As demonstrated in Yang, Mo & van den Bosch 
(2003) and van den Bosch, Yang & Mo (2003), the CLF is 
well constrained by the galaxy luminosity function and by 
the galaxy-galaxy correlation lengths as function of lumi- 
nosity. In the MGRSs used here we use the CLF with ID # 
6 given in Table 1 of van den Bosch et al. (2005). We have 
tested that none of our results depend significantly on this 
particular choice for the CLF. 

Because of the mass resolution of the simulations and 
because of the completeness limit of the 2dFGRS, we adopt 
a minimum galaxy luminosity of L m i n = 10 7 h~ 2 L©. The 
mean number of galaxies with L > L m i n that resides in a 
halo of mass M is given by 



(N) t 



$(L|M) dL 



(16) 



In order to Monte-Carlo sample occupation numbers for in- 
dividual haloes, one requires the full probability distribu- 
tion P(N\M) (with N an integer) of which (N)m gives the 
mean. We differentiate between satellite galaxies and cen- 
tral galaxies. The total number of galaxies per halo is the 
sum of iV C en , the number of central galaxies which is either 
one or zero, and iV S at, the (unlimited) number of satellite 
galaxies. We assume that iV sa t follows a Poisson distribu- 
tion and require that iV sa t = whenever A ccn = 0. The 
halo occupation distribution is thus specified as follows: if 
{N)m < 1 then A sat = and A cen is either zero (with prob- 
ability P = 1 — (N)m) or one (with probability P = (N)m)- 
If {N)m > 1 then A con = 1 and A sat is drawn from a Poisson 
distribution with a mean of {N)m — 1. 

We follow Yang et al. (2004) and draw the luminosity 
of the brightest galaxy in each halo from 3>(L|M) using the 
restriction that L > L\ with Li defined by 



$(L|M)dL = 1 . 



(17) 



The luminosities of the satellite galaxies are also drawn from 
$(L|M), but with the restriction L m i n < L < L\. 

Next we assign all galaxies a position and velocity 
within their halo, using the number density distributions 
and (isotropic) velocity dispersion profiles given in Section 4. 
Note that this implicitly assumes that all haloes, as well as 
their galaxy populations, are relaxed. Halo concentrations 
as function of halo mass are computed using the relation 
given by Eke, Navarro & Steinmetz (2001). 
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Table 1. Comparison between MGRSs and 2dFGRS 



MGRS 


b ve \ 


^rad 




/ccn 




Pks 150, 10001 


Pks 150, 2001 

rVO L * J 


Pks 1200, 3501 


Pks [350, 10001 


(1) 


(2) 


(3) 




(4) 




(5) 


(6) 


(7) 


(8) 


M .o 


0.0 


0.0 




0.0 




1.5 x 10~ 6 


2.6 x lO" 1 


5.1 x 10~ 4 


2.9 x 10~ 2 


M .i 


0.1 


2.8 x 10" 


-3 


7.3 x 10" 


-4 


6.0 x 10~ 4 


8.7 x 10" 1 


2.3 x lO" 1 


7.4 x 10~ 2 


M . 2 


0.2 


1.0 x 10- 


-2 


3.3 x 10" 


-3 


3.7 x 10~ 3 


4.9 x lO" 1 


6.2 x lO" 1 


8.7 x 10~ 2 


M . 3 


0.3 


2.4 x 10~ 


-2 


8.5 x 10" 


-3 


4.3 x 10~ 3 


4.3 x lO" 1 


2.4 x lO" 1 


1.6 x 10" 1 


M .4 


0.4 


4.2 x 10" 


-2 


1.8 x 10" 


-2 


1.8 x lO" 1 


8.2 x 10~ 2 


1.3 x lO" 1 


3.5 x lO" 1 


M . 5 


0.5 


7.2 x 10- 




3.5 x 10" 


-2 


2.8 x lO" 1 


7.5 x 10~ 2 


1.0 x 10- 1 


3.8 x 10- 1 


Mo. 6 


0.6 


1.1 x 10" 


-1 


6.6 x 10" 


-2 


8.5 x 10~ 3 


4.0 x 10~ 2 


1.1 x 10~ 3 


1.2 x 10" 1 


M . 7 


0.7 


1.8 x 10" 


-1 


1.3 x 10" 


-1 


3.3 x 10~ 3 


3.3 x 10~ 3 


1.2 x 10~ 2 


7.0 x 10~ 2 


M . 8 


0.8 


3.4 x 10~ 


■1 


3.3 x 10" 


-1 


1.0 x 10~ 5 


2.4 x 10~ 5 


1.1 x 10~ 3 


7.4 x 10~ 3 


Mi.o 


1.0 


1.0 








1.6 x 10~ 9 


1.9 x 10~ 5 


1.3 x 10~ 6 


5.3 x 10~ 6 



The MGRSs used for comparison with the 2dFGRS. Column (1) lists the ID of the MGRS. Columns (2), (3), and (4) list the velocity 
bias, spatial bias, and value of / CO n, respectively (see Section 4 for definitions). Finally, columns (5)-(8) list the KS probabilities Pks that 
the distributions of TZ extracted from these MGRS are consistent with those of the 2dFGRS for four different intervals in cr a , indicated 
by the values in square brackets (in kms" 1 ). 



5.2 Creating Mock Surveys 

The 2HFGRS uses a multifibre spectrograph to obtain red- 
shifts. However, because of the physical size of the fibers, 
when two galaxies are closer than ~ 30 arcsec in projec- 
tion only one of them can be targeted. Furthermore, due 
to clustering, some areas on the sky contain more galaxies 
within a single two-degree field than the available number of 
fibers. By using a sophisticated tiling strategy these prob- 
lems are largely overcome, yielding a fairly uniform sam- 
pling rate. Nevertheless, some spatial non-uniformities re- 
main. In addition, fainter galaxies yield noisier spectra, and 
therefore less accurate redshifts. All these effects combined 
result in a redshift completeness which depends on both po- 
sition on the sky and on apparent magnitude. The 2dFGRS 
team has constructed maps that parameterize this position 
and magnitude dependent completeness (Colless et al. 2001; 
Norberg et al. 2002), and which facilitate a simulation of 
these effects in our MGRSs. However, as it turns out, the 
completeness depends also on the angular separation, 0, be- 
tween galaxy pairs (see Hawkins et al. 2003). This is largely 
due to the problem of fiber collisions, which has not been 
completely corrected for by the tiling strategy. Finally, Nor- 
berg et al. (2002) have shown that the parent catalogue of 
the 2dFGRS, the APM catalogue, is only 91% complete. As 
shown in van den Bosch et al. (2005), this incompleteness 
is, at least partially, due to image blending in the APM cat- 
alogue (see also Cole et al. 2001). Based on this information 
we mimic the various observational selection and complete- 
ness effects in the 2dFGRS using the following steps: 

(i) We define a (a, (5)-coordinate frame with respect to 
the virtual observer at the center of the stack of simulation 
boxes, and remove all galaxies that are not located in the 
areas equivalent to the NGP and SGP regions of the 2dF- 
GRS. 

(ii) For each galaxy we compute the apparent magnitude 
according to its luminosity and distance, to which we add a 
rms error of 0.15 mag. Since galaxies in the 2dFGRS were 
pruned by apparent magnitude before a K-correction was 
applied, we proceed as follows: We first apply a negative 
K-correction, then select galaxies according to the position- 
dependent magnitude limit (obtained using the apparent 



magnitude limit masks provided by the 2dFGRS team) , and 
finally K-correct the magnitudes back to their rest-frame bj- 
band. Throughout we use the type-dependent K-corrections 
given in Madgwick et al. (2002). 

(iii) For each galaxy we compute the redshift as 'seen' by 
the virtual observer. We take the observational velocity un- 
certainties into account by adding a random velocity drawn 
from a Gaussian distribution with dispersion 85 kms" 1 . 

(iv) To take account of the position- and magnitude- 
dependent completeness of the 2dFGRS, we randomly sam- 
ple each galaxy using the completeness masks provided by 
the 2dFGRS team. 

(v) To take account of the fiber-collision induced incom- 
pleteness, we compute the angular separations 9 between 
all galaxy pairs and remove galaxies based on a probability 
p{9), which we tune (by trial and error) so that we reproduce 
the pair-separation incompleteness quantified by Hawkins et 
al. (2003). 

(vi) To take account of the incompleteness in the APM 
catalogue due to image blending we model the characteristic 
size of a galaxy as 

^ 1 = 15^ 1 k P c( I ^^ ?1 -) 1/3 (18) 

and define the critical projection angle # max = R SBi \/Da, 
with Da the angular diameter distance of the galaxy. We 
then remove the faintest galaxy from all pairs for which < 

^max • 

(vii) Finally, we remove a number of galaxies completely 
at random to bring the total fraction of removed galaxies, 
including those removed under (v) and (vi), to 9 percent. 

As shown in van den Bosch et al. (2005) , this procedure 
results in mock 2dFGRS catalogues that accurately mimic 
all the various incompleteness effects, allowing for a direct, 
one-to-one comparison with the true 2dFGRS. 



6 RESULTS 

Using the method outlined above, we construct a set of ten 
MGRSs that only differ in the value of 6 vc i. Table 1 lists these 
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Figure 5. The KS-probability that the 1Z distribution obtained from the 2dFGRS groups is consistent with that obtained from our 
MGRSs, as function of fe vc i. Results are shown for four cr s -intervals, indicated in square brackets in each panel. 



mocks together with their corresponding values of & ra( j and 
/con, computed for a halo with a concentration parameter 
c = 10. In the case of Mi.o, we deviated somewhat from the 
procedure described in Section 4. Rather than giving the 
central galaxies a probability distribution (10), we simply 
treat the brightest galaxy as a satellite galaxy so that 6 vc i = 
1.0. Note that in this case / CO n is not defined. For each of our 
ten MGRSs, we construct a group catalogue as described in 
Section 3.1, using those mock galaxies that are in the redshift 
range 0.01 < z < 0.20 and with a completeness c > 0.8 (this 
mimics our selection from the 2dFGRS). In what follows we 
restrict our analysis to groups with four or more members 
and with 50kms _1 < a s < 1000 kms -1 . 

To illustrate the importance of using MGRSs, the left- 
hand panel of Fig. 3 shows the cumulative distributions of 
\TZ\ obtained from the groups in Mo.o for which fe vc i = 0.0. 
In this mock all brightest halo galaxies have been located 
at rest at the center of the halo. The gray area indicates 
the area bounded by Student t-distributions with 3 and 9 



degrees of freedom (corresponding to systems with 4 and 
10 satellites, respectively, which spans the range covered 
by the vast majority of our groups). In principle, since this 
MGRS obeys the null-hypothesis of the CGP, the resulting 
P(< [R-\) should fall in this range. Clearly it doesn't, espe- 
cially not for groups with 50kms _1 < a s < 200kms~ 1 (solid 
line). This owes to the completeness effects in the survey, 
and to the fact that our group finder is not perfect and (un- 
avoidably) selects interlopers. As we discussed in Section 2, 
these effects systematically broaden P{TV) (see discussion 
in Section 2). Since the impact of one or two interlopers is 
much stronger in low mass groups, which have fewer mem- 
bers, the P(1Z) of groups with low a s deviates more from 
the predicted Student t-distribution than that of more mas- 
sive groups. This clearly demonstrates that one needs to 
take interlopers and completeness effects into account, in a 
statistical sense, when interpreting the distribution of 1Z ob- 
tained from the 2dFGRS. The MGRSs used here are ideally 
suited for this task. 
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The right-hand panel of Fig. 3 shows the same results 
as in the left-hand panel, but now based on Mi.o for which 
6 vo i = 1.0. Clearly, for this MGRS the P(7V) are significantly 
broader than for Mo.o- This demonstrates that, despite the 
interloper/completeness problem, the detailed distributions 
of TZ obtained from group catalogues do contain useful in- 
formation that we can use to constrain 6 ve i. 

In Fig. 4 we compare P(< \1Z\) obtained from our 2dF- 
GRS group catalogue (gray dots), to those obtained from 
three MGRSs with different values of 6 ve i, as indicated. The 
upper-left panel plots the results using all groups in the full 
range of a s considered. Clearly, the MGRS with & vo i = 0.0 
(i.e. the one that fulfills the null-hypothesis of the CGP) pre- 
dicts a narrower distribution of \1Z\ than that found for the 
2dFGRS. Although a model in which there is absolutely no 
segregation of the brightest galaxy, i.e., 6 ve i = 1.0, predicts a 
distribution that is clearly too broad, the intermediate case, 
with fe vc i = 0.5 matches the 2dFGRS results nicely. Similar, 
though somewhat more noisy results (because of the smaller 
number of groups involved) are obtained for the three sepa- 
rate bins of a s shown in the other three panels. 

Fig. 5 plots the KS probability, Pks, that the PQR\) of 
the 2dFGRS and the MGRS are drawn from the same dis- 
tribution, as function of 6 vo i (see also Table 1). This analy- 
sis shows that the brightest halo galaxies in the 2dFGRS 
(and thus also the SDSS) have a non-zero velocity with 
respect to the coordinate frame in which the mean satel- 
lite motion is zero. When using all groups in the range 
50 kms -1 < a s < 1000 kms -1 , the data is most consis- 
tent with a velocity bias of 6 vc i — 0.5, while the null hy- 
pothesis of the CGP is rejected at a high level of confidence 
(Pks = 1.5 x 10~ 6 , cf. Table 1). When analyzing the three 
a"s-subsamples, their is hint that the velocity bias of bright- 
est halo galaxies is more pronounced in more massive haloes. 
We caution, however, that these results are more noisy due 
to the smaller number statistics. 



7 DISCUSSION 

Our finding that brightest halo galaxies have, on average, a 
specific kinetic energy that is ~ 25 percent of that of satel- 
lite galaxies has two possible interpretations. First of all, 
the central galaxy may not be at rest with respect to the 
virialized dark matter halo. This scenario, which we here- 
after refer to as the Non-Relaxed Galaxy (NRG) scenario, is 
illustrated in the left-hand panel of Fig. 6. The right-hand 
panel depicts the second possible scenario; that of a Non- 
Relaxed Halo (NRH). In this case, the brightest halo galaxy 
is located at rest with respect to the minimum of the dark 
matter potential, but the dark matter mass distribution is 
not relaxed and reveals a clear m = 1 mode (i.e., the po- 
tential minimum does not coincide with the barycenter). In 
both scenarios, the brightest halo galaxy has a net velocity 
with respect to the coordinate frame defined by the mean 
motion of the satellites. Note that, although our MGRSs arc 
based on the NRG scenario, to first order it also mimics the 
NRH scenario, so that our comparison between MGRS and 
2dFGRS applies to both cases. 

Of the two scenarios illustrated in Fig. 6, the most likely 
one is the NRH scenario. It seems to fit naturally within 
a hierarchical picture of structure formation, where haloes 



continue to grow in mass by accretion and merging. It is 
also in accord with our finding that b vc i is larger in more 
massive haloes, which form later and are thus expected to 
be less relaxed. A potential problem for this scenario, how- 
ever, is the fact that in the ACDM concordance cosmology 
the growth rate of structures should drop fairly rapidly at 
the current epoch. Numerical simulations are ideally suited 
to investigate whether a value of b vc \ ~ 0.5 is a natural 
outcome of structure formation in the ACDM concordance 
cosmology or not. Using a combination of numerical simu- 
lations and semi-analytical models of galaxy formation in 
a ACDM cosmology, Diaferio et al. (1999) found that the 
central galaxy has an average velocity with respect to the 
halo barycenter of ~ 80kms~ . With a median halo mass of 
~ 10 13 ft _1 Mq, this corresponds to an average velocity bias 
of 6 vc i — 0.35, in reasonable agreement with our results. 
Yoshikawa, Jing & Borner (2003) used a smoothed-particle 
hydrodynamics (SPH) simulation of galaxy formation in a 
ACDM universe, and found that the average velocity differ- 
ence between the most massive halo galaxy and the mass 
center of the dark matter halo is (v c — Wd m ) ~ 0.5(<7dm) (see 
their Fig. 10) §. This is in excellent agreement with our best- 
fit value of & vc i =0.5. These two simulations suggest that our 
results are in perfect accord with the ACDM concordance 
cosmology. Note, however, that Yoshikawa et al. (2003) did 
not investigate whether the most massive halo galaxy has a 
net velocity with respect to the most bound halo particle, 
while in Diaferio et al. (1999) the central galaxy is associ- 
ated with the most bound halo particles by construction. 
Therefore, we not use either of these results to discriminate 
between the NRG and NRH scenarios. 

The NRG scenario appears unlikely at first sight, as dy- 
namical friction against the highly concentrated dark mat- 
ter halo should quickly damp any oscillatory motion. On 
the other hand, if dark matter haloes are cored, rather than 
cusped, the oscillations may persist for a much longer time 
(cf. Bontekoe 1988). This possibility is interesting in light 
of various independent claims for cored dark matter haloes, 
based on rotation curves of dwarf and low surface bright- 
ness galaxies (e.g., Moore 1994; Flores & Primack 1994; 
Borriello & Salucci 2001; de Blok et al. 2001; de Blok & 
Bosma 2002, but see also van den Bosch et al. 1999; van den 
Bosch & Swaters 2001; Dutton et al. 2005), on the observed 
pattern speeds of barred galaxies (Debattista & Sellwood 
1998, 2000), and on the longevity of the lopsidedness of disk 
galaxies (Levine & Sparke 1998). A more in-depth study of 
the damping rate of these kind of oscillations in dark mat- 
ter haloes, both cusped and cored, could shed more light on 
these issues. 

Independent of which of the aforementioned scenarios 
is responsible for the non-zero velocity bias of the brightest 
halo galaxies, it has important implications for various ar- 
eas in astrophysics. First of all, it has an important impact 
on the use of satellite kinematics to infer halo masses. Since 
the number of detectable satellites in individual systems is 
generally small, one typically stacks the data on many host- 
satellite pairs to obtain statistical estimates of halo masses 

§ We assumed isotropy to convert the three-dimensional velocity 
dispersion of the dark matter particles quoted in Yoshikawa et 
al. (2003) to the one-dimensional velocity dispersion (o"<j m ). 
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Figure 6. An illustration of the two different configurations that are both consistent with our inferred offset between the brightest 
halo galaxy and the satellite galaxies. Contours depict equipotentials of the dark matter haloes, while filled ellipses and circles represent 
brightest halo galaxies and satellite galaxies, respectively. In the Non-Relaxed Galaxy (NRG) scenario, the brightest halo galaxy oscillates 
in a fully relaxed halo. In the Non-Relaxed Halo (NRH) scenario, on the other hand, the central galaxy coincides with the minimum of 
the halo potential, but the centers of different equipotcntial surfaces are offset from each other. See text for a detailed discussion. 



(Erickson, Gottesman & Hunter 1987; Zaritsky et al. 1993, 
1997; Zaritsky & White 1994; McKay et al. 2002; Brainerd & 
Specian 2003; Prada et al. 2003; van den Bosch et al. 2004). 
The halo mass is typically derived from the dispersion, a cs , 
of the distribution of the velocity difference between host 
and satellite galaxies. This derivation rests on the (standard) 
assumptions that the host galaxies (i.e., the brightest halo 
galaxies) are at rest with respect to the center of a relaxed 
dark matter halo. If the satellite galaxies have the same kine- 
matics as dark matter particles, then ct cs = f"dm oc M 1 ^ 3 . 
However, in the case of the NRH scenario, one simply can 
not use (satellite) kinematics to infer reliable halo masses, 
as the crucial assumption of a virialized system is not cor- 
rect. In the case of the NRG scenario, on the other hand, 
the system is relaxed but, because of the non-zero velocity 
bias, we have that er cs = y/1 + ftveiCdm- If one were not to 
correct for 6 ve l, the inferred halo mass will be overestimated 
by a factor (l + 6 vc i) 3//2 (corresponding to ~ 1.85 in the case 
of 6 ve i = 0.5). 

The results presented here also have potentially impor- 
tant implications for (strong) gravitational lensing. In both 
the NRG and NRH scenarios one expects a strong, 'external' 
shear due to the dark matter halo, which should leave sig- 
natures in the image configurations and time delays of the 
lens. In fact, this 'external' shear may already have been 
detected. As shown in Keeton, Kochanek & Seljak (1997), 
fitting four-image lenses almost always requires an indepen- 
dent external shear that is not aligned with the light of the 
lens. Although this may reflect a misalignement between the 
luminous galaxy and dark matter halo, in agreement with 
the results presented here, there are alternative sources of 
external shear (nearby galaxies, large-scale structure along 
the line-of-sight) that may leave a similar signal in the lens 



configuration. A more thorough, systematic study of multi- 
ply lensed systems is therefore required to put constraints 
on the spatial bias of brightest halo galaxies. In fact, strong 
gravitational lensing is probably the only method that can 
be used to detect an offset between halo and galaxy in indi- 
vidual systems, and to discriminate between the NRH and 
NRG scenarios. 

A non-zero (r ccn ) also impacts on the internal structure 
and dynamics of central galaxies. As the galaxy oscillates in 
the dark matter halo (NRG scenario), or the halo relaxes 
around the central galaxy (NRH scenario), it is constantly 
subjected to tidal forces that may trigger bar-instabilities in 
otherwise stable disks, may cause excessive heating of the 
disk, and may create lopsidedness (Levine & Sparke 1998; 
Noordermeer, Sparke & Levine 2001). Detailed studies have 
revealed lopsidedness (either in the kinematics or the pho- 
tometry) in about half of all disk galaxies studied (e.g., 
Richter & Sancisi 1994; Zaritsky & Rix 1997; Rudnick & 
Rix 1998; Haynes et al. 1998; Matthews, van Driel & Gal- 
lagher 1998; Swaters 1999). In fact, as shown by Bissantz, 
Englmaier & Gerhard (2003), the morphology and kinemat- 
ics of gas in the inner few kpc of the Milky Way (in particular 
the 3kpc-arm) may indicate the presence of a similar m = 1 
asymmetry in our own galaxy (cf. Fux 1999). The high fre- 
quency of lopsided and barred disk galaxies therefore seems 
to be in support of a non-zero (Veen }, whether it reflects a 
non-relaxed halo or a non-relaxed galaxy. Taking our results 
at face value, it is clear that any study of disk stability that 
ignores these strong distortions and time-variability of the 
potential may be missing an essential ingredient. 

Finally, as mentioned in the introduction, a non-zero 
& V ei also plays a role in halo occupation models. Using the 
method described in detail in Yang et al. (2004), we com- 
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puted the projected two-point correlation function and pair- 
wise peculiar velocity dispersions of MGRSs Mo.o, Mo. 5, and 
M1.0. The differences are found to be extremely small, well 
below the errors due to cosmic variance. Therefore, for all 
practical purposes, it suffices to model the phase-space pa- 
rameters of galaxies in dark matter haloes with 6 vo i = (as 
is generally done), when computing galaxy-galaxy correla- 
tion functions based on halo occupation distributions. 



8 CONCLUSIONS 

According to the standard paradigm of structure formation, 
the brightest galaxy in a dark matter halo should reside at 
rest at the center of the potential well. In order to test this 
'central galaxy paradigm' (CGP), we used the halo-based 
galaxy group finder of Yang et al. (2005a) to construct group 
catalogues from the 2dFGRS and SDSS. For each group we 
compute the statistic 1Z, defined as the difference between 
the velocity of the brightest group galaxy and the average 
velocity of the other group members (satellites), normalized 
by the unbiased estimator of the velocity dispersion of the 
satellite galaxies. If the null-hypothesis of the CGP is cor- 
rect, 1Z should follow a Student ^distribution. If, on the 
other hand, brightest halo galaxies have a non-zero velocity 
bias with respect to the satellite galaxies, the ^-distribution 
should be significantly broader. The applicability of this '7Z- 
test' depends critically on how well one can group those 
galaxies that belong to the same dark matter halo. Although 
our group finder is well tested and calibrated, it is not per- 
fect, and unavoidably selects interloper galaxies as group 
members. In addition, redshift surveys suffer from various 
incompleteness effects. We have shown that these effects re- 
sult in a broadening of the ^-distribution, which, when not 
accounted for, may give the false impression that the CGP 
is ruled out. 

In order to take interlopers and incompleteness effects 
into account, and thus allow for a fair comparison with the 
data, we construct detailed mock galaxy redshift surveys 
that can be compared with the 2dFGRS on a one-to-one ba- 
sis. We apply our 7?.-statistic to the galaxy groups selected 
from these MGRS, which we compare to those obtained 
from the 2dFGRS using the Kolmogorov-Smirnov test. This 
shows that the CGP is inconsistent with the data at high 
confidence, and that instead the brightest halo galaxies have 
a specific kinetic energy that is about 25 percent of that of 
satellite galaxies. For a typical, relaxed, cold dark matter 
halo this corresponds to an expectation value for the offset 
between galaxy and halo of ~ 3 percent of the virial radius, 
comparable to the characteristic radius of the galaxy. In ad- 
dition, we find a weak hint that the velocity bias of brightest 
halo galaxies is larger in more massive haloes. 

We have focussed mainly on the ^-distributions ob- 
tained from the 2dFGRS, simply because we have accurate 
MGRSs available for this data set. However, we have shown 
that the ^-distributions obtained from groups in the SDSS 
are in excellent agreement with those obtained from the 2dF- 
GRS, suggesting that the SDSS is also inconsistent with the 
CGP. 

Undoubtedly, the most important implication of our re- 
sults is the puzzling question as to the origin of the offset 
between the central galaxy and its satellites. Probably the 



most likely explanation is that the majority of dark matter 
haloes are not yet fully relaxed. In this case, the brightest 
halo galaxy may still coincide with the minimum of the po- 
tential well, but that minimum does not coincide with the 
center of mass measured over the entire halo. Although this 
picture seems consistent with our finding that the specific 
kinetic energy of the brightest halo galaxies is larger in more 
massive haloes, which form later, detailed numerical simula- 
tions are required to investigate whether the typical growth 
rates of dark matter haloes are sufficiently large and violent 
to explain our findings in a ACDM concordance cosmology. 
A first hint that this is indeed the case comes from the simu- 
lations of Diaferio et al. (1999) and Yoshikawa et al. (2003), 
which reveal a velocity bias of brightest halo galaxies that is 
very similar to that found here. Although this suggests that 
a non-zero velocity bias is a natural outcome of structure 
formation in a ACDM cosmology, it still needs to be verified 
whether, in these simulations, the central galaxy is at rest 
with respect to the minimum of the potential well. 

An alternative explanation for the non-zero velocity of 
the brightest halo galaxies with respect to the satellite galax- 
ies may be that the halo is relaxed, but that the brightest 
halo galaxy oscillates in the central potential well. If the 
dark matter halo is strongly concentrated, as expected for 
typical cold dark matter haloes, one would naively expect 
that any such oscillation is quickly damped by dynamical 
friction. However, this damping timescale may be signifi- 
cantly longer if there is less dark matter in the center of the 
halo than anticipated; i.e., the density distribution is cored 
rather than cusped. This possibility is interesting in light 
of the recent claims for cored haloes based on the observed 
rotation curves and bar pattern speeds of disk galaxies. 

In either case, the brightest halo galaxy is expected to 
experience a time- varying tidal field. This strongly questions 
the applicability of (numerical) studies of galaxy dynamics, 
and in particular of stability analyzes that make the assump- 
tion that the galaxy is at rest at the center of a relaxed dark 
matter halo. In particular, it may explain the high frequency 
and longevity of bars and lopsidedness in disk galaxies. The 
fact that we find evidence for a non-zero velocity of the 
brightest halo galaxy with respect to the satellite galaxies 
also has important implications for the determination of halo 
masses based on the kinematics of host-satellite systems, and 
for the modeling of strong gravitational lenses. For the pur- 
pose of computing galaxy-galaxy correlation functions based 
on halo occupation models, however, one can safely ignore 
the fact that the CGP does not hold, and make the simple 
ansatz that the brightest halo galaxy resides at rest at the 
halo center. 
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